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James et al. presented simulations that appar- 
ently falsify the analytical result by Bastolla et al. 2 , 
who showed that nested mutualistic interactions de- 
crease interspecific competition and increase biodi- 
versity in model ecosystems. This contradiction, 
however, mainly stems from the incorrect applica- 
tion of formulas derived for fully connected net- 
works to empirical, sparse networks. 

Bastolla et al. 2 showed analytically that a model 
of mutualistic networks has two solutions, the weak 
regime and the strong regime. The former leads to 
a stable community up to a given threshold in the 
strength of mutualistic interactions. Beyond this 
threshold, the weak solution becomes unstable and 
the strong regime becomes stable. 

James et al. used what they thought to be a set 
of parameters within the weak regime. However, 
this set of parameters was determined by using the 
formula for a fully connected network. For empir- 
ical networks, the condition of the weak regime is 
that the effective competition matrix is positive def- 
inite, which is the case in only 22 out of the 59 
studied networks. Thus, the simulations of James 
et al. 1 occur mainly in another regime than the one 
assumed by the theoretical development of Bastolla 
et al. 2 For example, the notion of effective compe- 
tition is not even defined in the strong regime, al- 
ready casting doubts on the appropriateness of their 
comparison. 
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The same confusion of formulas that stand only 
for fully connected networks lead James et al. 1 not 
to include the effect of the productivity vector vari- 
ability on biodiversity. Consequently, these authors 
miscalculated the upper bound of biodiversity and 
missed the effects of nestedness 

Indeed, the general formula for the maximum bio- 
diversity (here extended for asymmetric effective 
competition matrices to mimic James et al.'s pa- 
rameter sampling procedure; see Methods) is not 
only a function of the effective competition, but also 
of the variability of the effective productivity vec- 
tor. This effective productivities are similar across 
species in fully connected networks, since they have 
the same number of interactions. Real networks, 
however, show a strong heterogeneity in the num- 
ber of interactions across species, and therefore in 
the components of their productivity vectors. Using 
the appropriate mathematical formulas (see Meth- 
ods), we confirm that the predicted upper bound 
of biodiversity is much closer to the observed val- 
ues. Specifically, the predicted value is 2 to 10 times 
larger than the observed value as opposed to 1000 
times larger as computed by James et al. (see Fig- 
ure la). 

Regarding the effect of nestedness on biodi- 
versity, our calculation with the extended cor- 
rect formula (see Methods) shows that the ef- 
fective competition — which is inversely related to 
biodiversity — is almost always negatively corre- 
lated to nestedness (Figure lb). Similarly, using 
a multilinear analysis (see Methods) we explicitly 
disentangle the effects of nestedness and those of 
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Figure 1: Ratio between the upper bound and actual 
biodiversity (a) and correlation between the effective 
competition and nestedness (b) using the generalized an- 
alytical expressions for structured networks. Parame- 
ter values as in Ref. 1 applied to the 22 mutualistic 
networks that are in the weak regime. Growth rates 
are sampled uniformly in [0.8, 1.4] to guarantee that the 
variability in the effective productivity vector is not too 
strongly correlated to nestedness. 1000 replicates are 
used. 



the effective productivity vector. We have done this 
independently for each network with the goal of re- 
moving the confounding effects of network size and 
connectance. This results in nestedness being al- 
ways significantly and positively related to biodi- 
versity. 

In summary, the apparently contradictory re- 
sults by James et al. can be turned upside-down 
if one applies the correct formulas for structured 
networks. Their approach would be like trying 
to prove that Pythagoras' theorem is false for the 
counter example of a non-right triangle. Both the 
expanded analytical results and the correct analysis 
of the simulations here shown confirm that the ob- 
served architecture of mutualistic networks indeed 
increases biodiversity. 

Methods 

We generalize the effective competition formula 
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diversity is generalized as 
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f?i (f9 2 ) is the angle between the effective produc- 
tivity vectors and the first left-(right)-singular vec- 
tor of the effective competition matrix. We relate 
biodiversity with nestedness and productivity as: 
log(biodiversity) ~ log (77) + log(l/A). 
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